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Abstract 

Intergalactic magnetic fields are assumed to have been spontaneously 
generated at the reheating stage of the early Universe, due to vacuum po- 
larization of non-Abelian gauge fields at high temperature. The fact that 
the screening mass of this type of fields has zero value was discovered re- 
cently. A procedure to estimate their field strengths, B{T), at different 
temperatures is here developed, and the value B{Tew) ~ lO^^G at the 
electroweak phase transition temperature, is derived by taking into con- 
sideration the present value of the intergalactic magnetic field strength, 
Bq ~ 10~^^G, coherent on the ~ 1 Mpc scale. As a particular case, the 
standard model is considered and the field scale at high temperature is es- 
timated in this case. Model dependent properties of the phenomena under 
investigation are briefly discussed, too. 



1 Introduction 

The recent experimental discovery of intergalactic magnetic fields having a field 
strength of the order B ~ lO^^^G is one of the most interesting events of 
modern cosmology [HE]- In Ref. [3] a model-independent, 95% CL interval 
1 X 10~^^G < B < 3 X 10~^^G was determined, and the femtoGauss values are 
actual field strengths in extragalactic space. This means, to start, that magnetic 
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fields are actually present everywhere in the Universe and influence various pro- 
cesses. Secondly, this renders most likely the cosmological origin of primordial 
magnetic flelds. From the theoretical viewpoint, this discovery restricts in an es- 
sential way the possible processes which may result in the creation of flelds in the 
hot Universe [2]. As a consequence, the search for mechanisms of fleld generation 
has intensifled. The most obvious candidates are primordial fluctuations, but 
there are other (for a review, see |l]-[6]). The challenge is to produce coherent 
magnetic flelds on very long scales in an almost empty intergalactic space. In this 
paper we will discuss some mechanism based on non-Abelian magnetic flelds. 

As was shown recently, spontaneous vacuum magnetization appears in non- 
Abelian gauge theories at high temperature. This was found by analytic methods 
in [7]- [To] and was conflrmed by means of lattice simulations in [11]. The basic 
idea rests on the known fact that spontaneous vacuum magnetization is the con- 
sequence of the spectrum of a color charged gluon, 

pl=pl + {2n + l)gB (n = -1, 0, 1, ...), (1) 

in a homogeneous magnetic background, B, described by the potential 

= Bx2 5^35^^ (2) 

where a is the weak isotopic index, and p\\ a momentum component along the 
fleld direction. Here, a tachyon mode is present in the ground state (n = — 1). 
In fact, one observes that J9q < 0, resulting from the interaction of the magnetic 
moment of the spin-1 charged particles with the magnetic fleld. This phenomenon 
was flrst discovered by Savvidy [12] at zero temperature, T = 0, and got known 
as the Savvidy vacuum. However, at zero and low temperature, this state is not 
stable. It decays under emission of gluons until the magnetic fleld B disappears. 
This picture changes with increasing temperature, where a stabilization sets in. 
The stabilization is due to vacuum polarization and depends on two dynamical 
parameters appearing for T ^ 0. These is a magnetic mass of the color charged 
gluon, nimagn., and an Ao-condensate, which is proportional to the Polyakov loop 
jl3] . This fleld conflguration is stable, its energy being below the perturbative 
one, and its minimum is reached for the fleld being of order gB ~ g^T'^/logT. 
This phenomenon is common to different SU{N) gauge flelds, which can be used 
to extend the standard SU{2) x U{l))ew x SU{3)c model of elementary particles. 

An important property of such temperature-dependent magnetic flelds is the 
vanishing of their magnetic mass, mmagn. = 0. This was found both in one-loop 
analytic calculations [H] and in lattice simulations [T5]. The mass parameter 
describes the inverse spatial scales of the transverse fleld components, similarly to 
the Debye mass mn, related to the inverse space scale for the electric (Coulomb) 
component. The absence of the screening mass means that the spontaneously 
generated Abelian chromomagnetic flelds are long range at high temperature, as 
is common for the U{1) magnetic fleld. Hence, it is reasonable to believe that. 
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at each stage of the evolution of the hot Universe, spontaneously created, strong, 
long-range magnetic fields of different types have been present. Since they are 
constant fields, their scales are coinciding with the horizon scale at each particular 
temperature. These fields may have critically influenced various processes and 
phase transitions. 

The idea of the present investigation is to relate the generation of intergalactic 
magnetic fields with the reheating epoch which immediately followed the infla- 
tionary stage of the universe evolution. At this epoch the temperature was the 
same in all volume of the expanded space (as it follows from the spectra of relic 
photons), and therefore, the field strength was also the same. Moreover, since 
both inflation and the vacuum polarization are causal processes, the coherence 
of the field at very large scales can be ensured. On the other hand, the creation 
of strong temperature dependent fields at this stage completely washed out the 
remnants of the fields generated at the inflation epoch. 

The dependence on the temperature of these fields differs from that of the 
usual U{1) magnetic fields. Recall that, in the latter case, as is commonly 
assumed the magnetic (hypermagnetic, in fact) field, created by some specific 
mechanism, is implemented in a hot plasma and evolves according to the law 
E ~ T^, which is a consequence of magnetic flux conservation (see, for instance, 
[6]). This scenario is actually considered, in magnetic hydrodynamics, when the 
evolution of large scale magnetic fields is investigated. The temperature depen- 
dence of a field of this type is ordinarily related with the high conductivity of the 
plasma. Thus, it is just a consequence of the classical theory of matter. How- 
ever, the classical theory is not always sufficient for describing the restored phase 
in models with spontaneous symmetry breaking. In this phase, quantum effects 
are essential. One of them is vacuum polarization, which results in spontaneous 
vacuum magnetization. Just this phenomenon regulates the dependence of the 
field strength on the temperature, and the magnetic flux is not conserved. In- 
stead, a specific flux value is generated at each temperature. This fact has to 
be taken into consideration when the cooling pattern of the hot plasma is inves- 
tigated at temperatures over that of the electroweak phase transition (EWPT), 
T > Tew ~ 100 GeV. This also concerns the SU{2)ew component of the elec- 
tromagnetic field. The classical theory is not able to account for it consistently. 
Thus, classical magnetic hydrodynamics could only start to work after the vac- 
uum magnetization has stopped. As we will show below, this does happen after 
the EWPT. 

It is worth to note that one of the difficult problems of magnetogenesis is to 
relate the field strengths generated via some given mechanism in the early Uni- 
verse with those of the present-day fields, what depends on numerous factors and 
is very model dependent. In magnetic hydrodynamics, the magnetic field evolu- 
tion was investigated in some detail for stochastic fields presumedly generated at 
either the electroweak or the confinement phase transitions (see [121 [IT] , and ref- 
erences therein). It was shown that their fate is strongly dependent on the scales 
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and initial conditions for these fields. In particular, magnetic flux conservation 
could not hold for some special conditions. We make a detailed comparison of 
these results with those coming from our approach, in what follows. An impor- 
tant point is that, in contrast with commonly investigated stochastic fields, we 
have here the possibility to deal with the evolution of a special solution of the 
field equations. Thus, its characteristics and properties are known and can be 
duely taken into consideration. 

In the frameworks of the standard model, we will here estimate the strength of 
the magnetic field at the temperature of the electroweak T^™ phase transition, as- 
suming that this field was spontaneously generated by a mechanism as described 
above. Although such phenomenon is nonperturbative, we carry out an actual 
calculation in the framework of a consistent effective potential (EP) accounting 
for the one-loop, V^'-^^ and the daisy (or ring), diagram contributions. In 

Sect. 2 we qualitatively describe, in quite more detail, the main aspects of the in- 
vestigated phenomena. In Sect. 3 the EP of an Abelian constant electromagnetic 
B field at finite temperature is obtained. It is used, in Sect. 4, to estimate the 
magnetic field strength at the EWPT temperatures. The scales of the magnetic 
fields at high temperatures, and their relation with the ones observed in inter- 
galactic space, are investigated in Sect. 5. A discussion of the results obtained, 
some conclusions, and prospects for further work are provided in the last section. 

2 Qualitative considerations 

In this section we describe, in a qualitative manner, the main aspects of the phe- 
nomena at issue. They are all consequences of asymptotic freedom and sponta- 
neous symmetry breaking at finite temperature. Our main assumption is that the 
intergalactic magnetic field has been spontaneously created at high temperature. 
We believe, this to be a fairly reasonable idea because, physically, magnetization 
is the consequence of a large magnetic moment for charged non-Abelian gauge 
fields (recall the gyromagnetic ratio 7 = 2 for 14^-bosons). This property results 
in the asymptotic freedom of the model in external fields. We discuss the proce- 
dure to relate the present value of the intergalactic magnetic field with the one 
generated in the restored phase. 

First, we note that, in non-Abelian gauge theories, magnetic fiux conservation 
does not hold at high temperatures, T > Tgw, where Tew is the temperature of 
EWPT. This is due to spontaneous vacuum magnetization, which depends on 
the temperature. The vacuum acts as a specific source generating classical fields. 
Second, the magnetization is strongly dependent on the scalar field condensate 
present in the vacuum at low temperature. This point was investigated at zero 
temperature by Goroku [IB]. For finite temperature, this is considered in the 
present paper for the first time. The observation is that, in both cases, the 
spontaneous vacuum magnetization takes place for small values of the scalar 
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field, 7^ 0, only. For the values of corresponding to any first order phase 
transition it does not happen. This means that, after the EWPT, the vacuum 
polarization ceases to generate magnetic fields, and magnetic flux conservation 
holds. As a result, the familiar dependence on the temperature, B ~ T^, is 
restored (for the fields spontaneously generated before the transition). However, 
that may be not the case, in general, for stochastic magnetic fields generated just 
at this transition [16]. One has to distinguish these types of fields, which have 
different nature and origin. 

Another aspect of the problem is the composite structure of the electromag- 
netic field Afj_. The potentials read 

where is the Z-boson potential, A^, are the Yang-Mills gauge field third 
projection in the weak isospin space and the potential of the hypercharge gauge 
fields, and g and g' are SU{2) and U{1)y couplings, correspondingly. After 
the electroweak phase transition, the Z-boson acquires a mass and the field is 
screened. Since the hypermagnetic field is not spontaneously generated, only the 
component = g' J^^j ^fg^^^T^ = sin6'w^^ is present at high temperature. 
Here 6'w is the Weinberg angle, tan6'w = g'/g- This is the only component 
responsible for the intergalactic magnetic field at low temperature. In the restored 
phase, b^ = 0, the complete weak-isospin chromomagnetic field A^^^ is unscreened. 
This is because the magnetic mass of this field is zero [15]. Thus, the field is of 
long range and this provides the coherence length Xb{T) to be sufficiently large. 
In fact, because the field is a constant, it has to cover all the horizon scale at 
the given temperature, Ab(T) ~ Rh(t)- This property is very important for our 
scenario of intergalactic magnetic field generation. Its cosmological consequences 
will be discussed below. In particular, magnetic fields of different types (color 
SU{?>), and others) can be spontaneously generated at high temperatures. 

Here, we continue with the description of the general field behavior related 
with the EWPT. In the restored phase, a scalar field condensate = 0, and the 
constituent of the weak isospin field corresponding to the magnetic one, is given 
by the expression 

i?(T) = sine^(T)i?(3)(T), (4) 

where B^^^{T) is the strength of the field generated spontaneously. After the 
phase transition, the scalar condensate is 7^ 0, and the field is partially screened. 

To relate the present value of the intergalactic magnetic field with the field 
which existed before the EWPT, we take into consideration that, after the phase 
transition, spontaneous vacuum magnetization does not take place. This property 
is derived in the next sections. Therefore, for the EWPT temperature Tew, we 
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can write 



Bo 



Bo 



(5) 



Here, Bq is the present value of the intergalactic magnetic field strength Bq ~ 
10~^^G. The left-hand-side relates the value B{Tew) with Bq. The right-hand- 
side gives a possibility to express the weak isospin magnetic field in the restored 
phase through Bq, knowing the temperature dependence of the Weinberg angle 
6w{T). This relation contains an arbitrary temperature normalization parameter 
T. It can be fixed for a given temperature and Bq. After that, the field strength 
values at various temperatures can be calculated. In particular, the total weak 
isospin field strength is given by the sum cos 9w{Teu,)B^^\Teu,) + B(Teu,)- 

Note also that the relation ([5]) is a consequence of the assumption that, for 
the field spontaneously generated at high temperature, the magnetic fiux con- 
servation holds, in the early Universe, after the EWPT. This means that the 
field is "frozen" in the plasma at large scales, and the magnetic turbulence pro- 
cesses do not affect that behavior. Although this is a most simple assumption, 
it requires a detailed discussion encompassing the results obtained recently in 
magnetic hydrodynamics (see Refs. fl6[ [T7] . and references therein). Details on 
the magneto-hydrodynamical processes occurring in the early universe can be 
found in numerous publications (see, for references, the review papers [SI E]). 
Here, we will just mention the main points which are important for our future 
considerations. The "frozen in" conditions are always realized for magnetic fields 
having the scales larger than the largest turbulence eddies. After a free decay 
stage of the magnetized plasma evolution, the field can be considered to be non- 
affected by turbulence [IT]. In connection with these results, it is clear that the 
magnetic fields generated at the EWPT are not sufficient, as such, to produce 
long-range correlated fields, and some additional processes must be included. 
This is because, even a field having the scale of the Hubble radius at temperature 
Tew ~ 100 GeV, is in fact correlated at the comoving scale Iq ~ 10~^ Mpc JT7\. 

But, in the case of a temperature-dependent vacuum polarization, other pos- 
sibilities exist. To obtain a large-scale correlated intergalactic magnetic field one 
can take into account that the temperature in the Universe is the same, to high 
accuracy, and of the order ST ~ 10~^. This is ensured by the reheating stage 
coming just after infiation. The temperatures at this stage are estimated to be of 
the order 10^^ — 10^^ GeV depending on the specific infiation model [H]. Hence, 
with high accuracy, the magnetic field strengths B(T) generated due to vacuum 
polarization are also the same in all regions of the universe. Since they are pro- 
duced by a causal process, as infiation is, a correlation of the fields at very large 
scale will result. In this way, coherent magnetic fields on huge scales could be 
generated. Such kind of fields have much larger scales than any turbulence eddy, 
therefore, they will not be infiuenced by turbulent processes. Thus, it is rea- 
sonable to assume that these magnetic fields were frozen into plasma after the 



EWPT. 
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An important aspect of the present scenario is that the knowledge of the 
particular theory, still unknown today, which may prove in future experiments 
to be the right extension to the standard model, is not that very important for 
estimating the field strength B at temperatures close to Tew This is so because 
the new gauge fields of the extended model will remain screened, at the high 
temperatures corresponding to the spontaneous symmetry breaking of some of 
the basic symmetries. At very high temperatures, when these symmetries will 
be restored, the corresponding magnetic fields do will exist. In this way the 
value of the field strength at the Planck era can be estimated (see Ref. |20], for 
comparison) . 



3 Effective potential at high temperature 

As we noted above, the spontaneous vacuum magnetization and the absence of 
the magnetic mass for the Abelian magnetic fields are nonperturbative effects 
to be determined, in particular, in lattice simulations [HI [15]. The main con- 
clusions of these investigations are that a stable magnetized vacuum does exist 
at high temperature and that the magnetic mass of the created field is zero. 
Concerning the actual value of the field strength, is is close to the one calcu- 
lated within the consistent effective potential which takes into account one-loop 
plus daisy diagrams. Thus, in the present investigation we restrict ourselves to 
such approximation. This is mainly with the purpose to be able to use analytic 
calculations in order to properly interpret the results. 

The complete EP for the standard model is given in the review [21]. In the 
present investigation we are interested in two limits, 

1. weak magnetic field and large scalar field condensate, h = eB/M^^ < 0^, 
= 0e/0o, /3 = 

2. the case of the restored sjTumetry, = 0, gB 7^ 0, T 7^ 0. 

For the first case, we show the absence of spontaneous vacuum magnetization 
at finite temperature. For the second, we estimate the field strength at high 
temperature. Here is the VT-boson mass at zero temperature, 0c a scalar 
field condensate, and 0o its value at zero temperature. 

To demonstrate the first property we consider the one-loop contribution of 
W^-bosons (see also Eq. (27) of Ref. p]). 



h °° 



n 

Here n labels discrete energy values and Ki{z) is the MacDonald function. 



(6) 
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The main goal of our investigation is the restored phase of the standard model. 
So, we obtain the high temperature contribution of the complete effective poten- 
tial relevant for this case using the results in Ref. [21]. First, we write down the 
one-loop ly-boson contribution as the sum of the pure Yang-Mills weak-isospin 
part {B = 5(3)), 

- ^^-^^^ + 0{9'B% (7) 
where r is a temperature normalization point, and the charged scalars [9], 

V}^\B,T) = -^4^^1og^ + i^^^^ + 0(/5^), (8) 



describing the contribution of longitudinal vector components. The first term in 
Eq. ([7j) is the tree-level energy of the field. This representation is convenient for 
the case of extended models including other gauge and scalar fields. Depending 
on the specific case, one can take into consideration the parts ([7j) and ([8]), corre- 
spondingly. In the standard model, the contribution of Eq. ^ has to be taken 
with a factor 2, due to the two charged scalar fields entering the scalar doublet 
of the model. In the case of the Two-Higgs-Doublet standard model, this factor 
must be 4, etc. The imaginary part is generated because of the unstable mode 
in the spectrum ([1]). It is canceled by the term appearing in the contribution of 
the daisy diagrams for the unstable mode [TO] . 



Vunstabie = ^[U{B,T,n = -1) - gB]'/' + (9) 

Here Il{B,T,n = —1) is the mean value for the charged gluon polarization ten- 
sor taken in the ground state n = — 1 of the spectrum ([1]). If this value is 
sufficiently large, spectrum stabilization due to radiation correction takes place. 
This possibility formally follows from the temperature and field dependence of 
the polarization tensor in the high temperature limit T — )■ oo [22]: 11(5, T,n = 

— 1) = c g^T\l gB, where c > is a constant which must be calculated explicitly. 
At high temperature the first term can be larger then gB. From Eqs. and 
([9]) it follows that the imaginary part is canceled. Hence, we see that having ac- 
counted for rings leads to vacuum stabilization even if H(i?, T,n = —1) is smaller 
then gB. Really, in the latter case, the imaginary part will be smaller than in 
Eq. ([7]). The high temperature limit of the fermion contribution looks as 

(y. 1 T 

Vfermion = Y] ^I'jB^ log (10) 

71 ^ T 
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where the sum is extended to all leptons and quarks, and g/ is the fermion 
electric charge in positron units. Hence, it follows that in the restored phase all 
the fermions give the same contribution. 

Let us now present the EP for ring diagrams describing the long range corre- 
lation corrections at finite temperature [231 [2^ , 

^ ~ )(0)-T7^Tr[noo(0)p/' 



^ing 



001 



+ 



24/32 

(noo(o))^ 

327r2 



127r/3 

47r 3 ■ 

°^/3(noo(0))V2 + 4~^ 



(11) 



where the trace means summation over all the contributing states, IIoo = n^(A; = 
0,T,B) for the Higgs particle; m|, = IIoo = ^oo{k = 0,T,B) are the zero- 
zero components of the polarization functions of gauge fields in the magnetic 
field taken at zero momenta, called the Debye mass squared, and 7 is Euler's 
gamma. These terms are of order ~ g^{^^^'^) in the coupling constants. The 
detailed calculation of these functions is given in Ref. [21]. We give the results 
for completeness, 

1 3e^ . 



n^(o) 



+ 



24/32 
2a 



(6A + 



sin2(2^^ 



sin 6„ 



(12) 



+ 



3/32 

Svrsin^ Oy^P^ 



\qfB\Kj 



^(3V2C(--,-)) 



Here Kf = 



^Yukawa 



and A is the scalar field coupling. The terms ~ T"^ 



yield standard contributions to the temperature mass squared coming from the 
boson and fermion sectors. The 5-dependent terms are negative (note the value 
of 3a/2C(— 1, |) = —0.39). They decrease the value of the screening mass at high 
temperature. The Debye masses squared for the photons, Z-bosons, and neutral 
current contributions are, respectively. 



m 



Da 



sin^ Q 



w 



3/3^ 



+ 0(eS/32 



( tan^ 6'w + 



1 



m 



g' 



D.neutr. 



4 cos2 6'w 
1 + 4 sin^ B 



3/3^ 



w 



+ 0(65/3^ 



+ 0{eBf 



8/32 cos2 

As one can see, the dependence on B appears at order 0{T^'^\ 
The PF-boson contribution to the Debye mass of the photons is 

- ^ ( 1 (^sin^w5)i/2\ 



(13) 



"ig' sin^ d 



w 



3/32 



27r/3 



(14) 
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An interesting feature of this expression is the negative sign of the next-to-leading 
terms which dependen on the field strength. Finally, we give the contribution of 
the high temperature part in Eq. ([9]) Il{B,T,n = —1) [21], 

U{B, T,n = -1) = Jasine^By/^ ^^^.as + 4z) . (15) 

This expression was calculated from the one-loop ly-boson polarization tensor 
in the external field at high temperature. It contains the imaginary part which 
comes from the unstable mode in the spectrum ([T]). Its value is small as compared 
to the real one. It is of the order of the usual damping constants in plasma at 
high temperature. It will be thus ignored in actual calculations in what follows. 
In fact, this part must be calculated in a more consistent scheme which starts 
with a regularized stable spectrum. On the other hand, as we noted above, the 
stability problem is a non-perturbative one. Stabilization can be realized not only 
through radiation corrections but also by some other mechanisms. For example, 
due to Aq condensation [7] at high temperature. A stable vacuum state was 
observed in lattice simulations [11], therefore, we believe that this problem has a 
positive solution. Summing up, we have all necessary ingredients to investigate 
the problem of interest. 



4 Magnetic field strength at 



ew 



Let us now show that the spontaneous vacuum magnetization does not take 
place at finite temperature and for non-small values of the scalar field condensate 
7^ 0. To this end we notice that the magnetization is produced by the gauge 
field contribution, given in Eq. ([6]). So, we consider the limit of gB/T"^ <^ 1 and 
0^ > h. For this case we use the asymptotic expansion of Ki{z), 

Mz)r.J^e-^(l + l--^ + ..), (16) 
\ 2z \ 8z 128z^ / 

where z = n(3{(l)'^ ± hY^"^. Now, we investigate the limit of /3 — )■ cx3, T/(f) <^ 1, 
where the leading contribution is given by the first term of the temperature sum 
in Eq. ([6]). We can also substitute (0^±/i)^/^ = 0(1 ± 2^). In this approximation, 
the sum of the tree level energy and (jS]) reads 



The second term is exponentially small and the stationary equation dV /dh = 
has the trivial solution h = 0. This estimate can be easily verified in a numerical 
calculation of the total effective potential. Hence, we conclude that after symme- 
try breaking the spontaneous vacuum magnetization does not take place, as was 
the case at zero temperature [H 
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To estimate the magnetic field strength in the restored phase at the EWPT 
temperature the total EP deduced in the previous section must be used and the 
parameters entering Eq. ([5]) need to be calculated. This can be best done numeri- 
cally. To explain the procedure, we consider here a part of this potential account- 
ing for the one-loop ly-boson contributions. The high temperature expansion for 
the EP coming from charged vector fields is given in Eq. ([7]) . Assuming stability 
of the vacuum state, we calculate the value of the chromomagnetic weak isospin 
field spontaneously generated at high temperature from Eqs. ([7]) and ([8]): 

1 o3 7-2 



This expression (and the complete one accounting for all contributions) gives the 
field strength at any temperature, T > Tg^. These formulas can be obtained for 
different types of particles. Before reducing to a specific value for it, we describe 
how to connect this expression with the intergalactic magnetic field Bq. We first 
relate the expression f|T8l) with an electromagnetic field after symmetry breaking, 
and then take into account the scales of the fields. 

Let us introduce the standard parameters and definitions, = g^/^Air), 
a = a„sin6'^, ay = {g')^/{4:n) and ta.ia^ 9y/(T) = ayiT) / a^{T), where a is the 
fine structure constant. To find the temperature dependence of the Weinberg 
angle, the behavior of the hypercharge coupling g' on the temperature has to 
be computed. From Eq. ([8]) it follows that this behavior is nontrivial. The 
logarithmic temperature-dependent term is negative. But, as is well known, in 
the asymptotically free models this sign must be changed into a positive value, 
due to the contributions of other fields. This particular value is model dependent. 
We will not calculate it in the present paper. Instead, for a rough estimate, we 
substitute the zero-temperature value: sin^6'w(r) = sin^6'w(0) = 0.23. 

For the given temperature of the EWPT, Te^, the magnetic field is 

S(Te^) = = sin0w(Te„)5(Te„). (19) 

Assuming T^^ = lOOGeV = W^eV and Tq = 2.7K = 2.3267 ■ 10~^eV, we obtain 

S(Tew) ~ 1.85 lO^^G. (20) 

This value can be considered as an estimate of the magnetic field strength at the 
EWPT. Hence, for the value of X = log we have the equation 

1 

^° = 2^^vw^ (TT^^^Vxp- ^^^^ 

Since all the values are known, logr can be computed. After that, the field 
strengths at different higher temperatures can be found. In fact, the main point 
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in obtaining these results is the assumption of magnetic flux conservation, as is 
frozen in a plasma. Information of a particular model is implemented in the factor 
sm9u;{Tew) in Eq. ( !T9l) . Needles to say, our estimate is a rough one, because we 
have ignored the temperature dependence of the Weinberg angle. To guess the 
value of the parameter r we take the field strength Bq ~ lO"^^, usually used in 
cosmology (see, e.g., [20]). In this case, from Eq. f l2T|) we obtain r ~ 300e\^. For 
the present-day value Bq ~ lO'^^G this parameter is much smaller. 

To take into account the fermion contribution Eq. ffTOl) . we have to substitute 
the expression j2^log'^ in Eq. ffTSl) and also in fl2T|) . with the value 

(5-i:i.|)>7- P^' 

In the above estimate, we have accounted for the one-loop part of the EP of order 
~ in the coupling constant. The ring diagrams have order ~ and give a 
small numeric correction to this result. As was mentioned before, had we taken 
into account all the terms listed in the previous section, the results would have 
not changed essentially. 

Let us compare now the value of the field strength fl20l) with the one calculated 
directly from the EP for the standard model in Ref. p5]. From Fig. 1 and Tab. 1 
of that paper, we find 

B^^'{T,^) ~ lO^OG, (23) 

what is much larger than the value (l20l) and just corresponds to the value of the 
present comoving field strength Bq ~ 10~^G. Note that this value was used in 
numerous investigations of magnetic fields in the early universe, before the recent 
discoveries [l]-[3]. Let us stress again that the field strength at higher temper- 
atures will depend on the particular model extending the standard one. Spon- 
taneous vacuum magnetization in the minimal supersymmetric standard model 
has been investigated in Ref. [21], and the field strength generated in this model 
is smaller as compared to the situation here considered. Also, Pollock [20] has 
investigated this problem for the case of the Planck era, where magnetic fields of 
the order B ~ lO^^G have been estimated. 

5 Magnetic field scale 

We now discuss in brief the scale of the field generated in the restored phase. 
This is a key point in relating expressions like Eqs. (fTS!) or (l23ll with Bq. In our 
consideration, the "frozen in" condition was used. Therefore, we are going to 
discuss its applicability in more detail. Note first that, if one assumes that after 
the EWPT the constant field B{Tew) was frozen in the plasma at the Hubble scale, 
Rh{T(,w), then its comoving coherence scale at present will be \b{Tq) = 6 ■ 10~^ 
pc [17]. This is much smaller than is needed. 
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We propose two, in fact related, ways to overcome such difficulty. The ffist is 
to take into consideration the reheating stage of the universe evolution. Accord- 
ing to the concepts of modern cosmology [19], this stage has existed just after 
inflation and is related with the latter causal stage. Just due to causality, the 
temperature in the universe after this stage is the same, in all domains of space, 
which could even be uncorrelated in later moments of time. Hence, at a given 
high temperature, T, the magnetic field generated due to vacuum polarization 
has the same strength B{T) everywhere in the universe. Formally, they could 
have different directions, in either external or internal spaces, although this point 
requires additional consideration. Different kind of chromomagnetic fields, of the 
type as in Eq. (|2]), can be spontaneously generated. Their nature depends on the 
particular model considered and is therefore unknown, as of now. But this is not 
essential for our consideration, here. The magnetic fields coherent on huge scales 
are expected to have been present in the early Universe. The origin of this coher- 
ence is ensured by the properties of the solution to the field equations (Eq. ([2])) 
and by causality at the inflationary epoch. The scales of the coherent field do- 
mains could be estimated on grounds of the gauge invariance. This idea, due to 
Feynman, was put in force in gluodynamics with the goal to determine possible 
magnetic vacuum structures |26]. Namely, to find a gauge invariant vacuum, on 
the basis of gauge non- invariant solutions (such as Eq. (EJ), one can consider a 
domain structure ensuring gauge invariance when a corresponding boundary is 
going around. This point requires further investigation. 

Most of the fields generated in the early Universe decouple and are screened 
at some energy (temperature) scales, when the corresponding scalar condensates 
have broken the background symmetries. So, the only unbroken symmetry at the 
EWPT remains the SU{2)eyj x f/(l) one. After the EWPT, when spontaneous 
magnetization stops, this field cannot be included in turbulent processes gener- 
ated by the transition. This is because the scale of the field, Eq. (|2]), is already 
much larger than any largest eddy of turbulence. As is usually believed, the 
size of a typical eddy is estimated as the inverse mass of the particles appearing 
after the transition [27]. Thus, the field evolves in accordance with the metric 
expansion and is implemented in a hot plasma, thus fulfilling the magnetic flux 
conservation law. And it finally results in the present day intergalactic mag- 
netic fields which could be correlated on ~ 1 Mpc scales. Note that an essential 
information on the processes that take place after the EWPT, obtained in the 
framework of magneto- hydrodynamics, is given in Refs. [T7] and [2S]. We will 
discuss the facts which have relevance to our problem in the last section. 

Another possible scenario is based on the stochastic processes considered al- 
ready by Hogan [29j in connection with the magnetic fields generated at ffist-order 
EWPT. A possible mechanism of field generation in that case was proposed by 
Vachaspati [27j. In the former paper, it was pointed out that magnetic fields 
correlated on large scales can be produced not only through causal processes 
but also by a stochastic random walk mechanism, if the magnetic lines gener- 
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ated in some domain of space "forget" about their origin. The field strength 
developed on large scales by this process (due to "straightening" of entangled 
magnetic fields) can be estimated as Bj^ ~ B/\/N, where N counts the number 
of domains, with the field i? of a given size, crossed by a magnetic line. The 
correlation length in this case can be much larger than the Rh{T). It can be 
estimated as Xb(T) ~ NRniT). In the paper [29], it has been also noticed that 
this mechanism is not applicable to the early universe, the reason being because 
magnetic lines do not penetrate freely though the plasma. This is really the case, 
if the properties of the plasma as such are taken into account. However, this is 
not the case if spontaneous vacuum magnetization occurs. Actually, at a given 
temperature, each uncorrelated domain of space having a Hubble radius Rh(T) 
is filled up with a constant magnetic field B(T), described by the potential (EJ. 
Its orientation in both external and internal spaces is arbitrary. Hence, a stochas- 
tic behavior of the field lines and the appearance of magnetic fields having large 
correlation lengths \b(T) > Rh(T) are expected. After the EWPT, these fields 
evolve as in the previous case. 

Note that, in both scenarios, all the fields generated at the inflation epoch are 
washed out by the vacuum polarization and leave no remnants at present. The 
reheating stage becomes more important. In closing this section, we also notice 
that the long range nature of the Abelian spontaneously generated magnetic fields 
is ensured by their zero magnetic mass (see [13 E]), what renders these fields 
unscreened, as is the case for usual U{1) fields. 

6 Discussion and conclusions 

We here summarize our main results. The key point in the problem under inves- 
tigation is the spontaneous vacuum magnetization process, which eliminates the 
magnetic flux conservation principle at high temperature. Vacuum polarization 
is responsible for the value of the field strength B{T) at each temperature and 
serves as its source. In fact, it converts heat into an ordered coherent state, which 
was not taken into consideration in previous studies of the early Universe. We 
have also shown here that, at finite temperature and after symmetry breaking, a 
scalar field condensate suppresses the magnetization. Hence, it follows that the 
actual nature of the model extending the standard one is not that essential at 
sufficiently low temperatures, when the decoupling of the heavy gauge fields has 
occurred. From this one can conclude, in particular, that the vacuum polarization 
"washes out" the relics of the magnetic fields generated at very high temperature 
or at the inflation stage. These statements are new and come as an interesting 
surprise, as compared with the standard notions based on the ubiquitous scenario 
with magnetic flux conservation. 

The present value of the intergalactic magnetic field strength has been here 
related with fields at high temperatures in the restored phase. Because of the 
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zero magnetic mass for Abelian magnetic fields (as discovered recently [HI [T5]). 
there is no problem for the generation of fields having large coherence scales. 
In our estimates, we have assumed that, basically, the field is of the order of 
the horizon scale, Xb(T) ~ Rh(t)- This seems reasonable because, at a given 
temperature, the field B{T) = const, generated due to vacuum polarization, 
occupies all space. Then, at the reheating stage (due to causality, present already 
at inflation), coherence on scales exceeding that of the horizon can be produced. 
Such large scale fields are not influenced by turbulent processes happening after 
the EWPT. They are frozen in the plasma and evolve according to the magnetic 
flux conservation law. In this scenario, a large scale domain structure is also 
permissible, what requires additional consideration. 

Knowing the particular properties of the extended model, it is possible to 
estimate the field strengths at any temperature. This can be done for different 
schemes of spontaneous symmetry breaking (restoration) by taking into account 
the fact that, after the decoupling of some massive gauge fields, the correspond- 
ing magnetic fields are screened. Thus, the higher the temperature, the larger 
the number of strong long-range magnetic fields of different types that will be 
generated in the early Universe. 

As we have found above, the field strengths at the EWPT temperature, esti- 
mated with account to the present-day value of the intergalactic magnetic field 
strength, ~ 10^^^ G, Eq. ( 120|) . or either directly from the vacuum magnetization 
in the standard model, Eq. ( 123|) . differ in six orders of magnitude. This huge 
deviation can be explained by the different scales of the fields considered. Let 
us check this possibility by using the second of the scenarios proposed in the 
previous section, for large scale field generation. Making use of the usual relation 
between the scale factor and the temperature, 

%^ = ^, (24) 

taken at the EWPT epoch, and the present-day parameters, T^w = 100 GeV = 
10^^ eV,To = 2.3267-10-^ eV. If one assumes that Xb{T} ^ a(T), then from (12^1). 
it follows that Xb{Tq) = 6 ■ 10~^ ps (see, for instance, [H]). On the other hand, if 
one takes Xb{To) = 1 Mpc, the value ^BiTew) = 2.33 ■ 10~^^ Mpc is obtained. At 
the same time, the horizon size is a(Te^) = 1.27 ■ 10~^^ Mpc, thus, XB{Tew) » 
a(Tew). Now, following an idea of Hogan P2], we relate the size of the correlated 
field with the random walk process. At T^^, we have XB{Tew) = Na{Tf,w), hence, 
we get roughly = 3-10^, and for the field strength "straightened" on the N- 
domain scale, ^ ^^^^ Therefore, accounting for the field strength value 
calculated for the standard model, Eq. ( |23l) . we obtain BisiTew) ~ 3 ■ 10^^ G (the 
subscript in Bis means "large scale"). This value is close the the value BisiTew) ~ 
2 ■ 10^^ G estimated in Eq. fl20|) . The remaining discrepancy can be explained in 
two ways. First, and obviously, as due to the roughness of our estimate. Second, 
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and more radically, by the necessity of substituting the standard model with 
another one. The latter point will be discussed in more detail below. 

Let us note a number of properties of the field under consideration, and com- 
pare them with the ones usually discussed in applications of magnetohydrody- 
namics to the early Universe. Here, we follow Refs. pT] and [28] which are close 
to our analysis. First, let us mention that the field energy density, pb = \, is 
proportional to g^T'^, what is much smaller than the radiation energy density, 
~ T^. Thus, the BBN condition (see [IZ]), Pb/ prad. « 1, is fulfilled. Second, as 
numerical simulations show , the turbulent process in the early Universe with 
magnetic field included is slower then in the laboratory. Turbulence can include 
a large scale field at the level of the largest eddies. For large scale fields, the 
free decay stage is important. At this epoch, which is strongly dependent on the 
initial conditions [28], turbulence is significantly decreasing, and after this very 
brief stage the field is not affected by turbulence any more. It is just frozen in the 
plasma. As we have shown, the spontaneous vacuum magnetization is stopped 
when the first-order phase transition ends. Fields of this type cannot be influ- 
enced by turbulence and it is thus reasonable to believe that after the EWPT 
the field evolves according to Eq. ([5]). Note also that these fields are non- helical 
ones. 

Our analysis has shown that, at the EWPT temperature, magnetic fields of 
the order B{T(.w) ~ lO^^G did very likely exist. To estimate their field strengths 
at higher temperatures, one has to take into consideration a number of features 
proper to the standard model and its particular extension at play. First, we note 
that quarks possess both electric and color charges. Therefore, there is a mixing 
between the color and usual magnetic fields owing to the quark loops. Second, 
there are peculiarities related with the particular content of the extended mod- 
els. For example, in the Two-Higgs-Doublet standard model the contribution 
~ {gBY^'^T in Eq. ([7]) is exactly canceled by the corresponding term in Eq. (|8]), 
because of the four charged scalar fields entering the model. They interact with 
gauge fields with the same coupling constant. However, in this model the dou- 
blets interact differently with fermions. This changes the effective couplings of 
the doublets with the gauge fields and results in non-complete cancelations. As 
a result, a suppression of the spontaneously created magnetic field is expected in 
this model. In principle, one should be able to explain, in this way, the discrep- 
ancy in the field strengths as discussed above. There can be other peculiarities 
which may influence the high temperature phase of the universe. They will re- 
quire further investigation and thus we leave this issue for a future publication. 
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